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Abstract 

A two dimensional model is introduced to study pattern formation, secondary 
instabilities and the transition to spatiotemporal chaos (weak turbulence) in 
parametric surface waves. The stability of a periodic standing wave state 
above onset is studied against Eckhaus, zig-zag and transverse amplitude 
modulations (TAM) as a function of the control parameter e and the detuning. 
A mechanism leading to a finite threshold for the TAM instability is identified. 
Numerical solutions of the model are in agreement with the stability diagram, 
and also reveal the existence of a transition to spatiotemporal chaotic states 
at a finite e. Power spectra of temporal fluctuations in the chaotic state are 
broadband, decaying as a power law of the frequency u}~^ with z w 4.0. 



PACS Numbers: 47.52.+j, 47.35.+i, 47.54.+r 
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A layer of incompressible fluid that is driven by a sinusoidal force normal to the free 
surface at rest exhibits parametrically excited surface waves, also known as Faraday waves 
||1|] . We introduce a two dimensional model that possesses many of the features observed in 
experimental studies of Faraday waves: a primary instability to a standing wave pattern, 
a secondary instability (transverse amplitude modulation or TAM) at larger amplitudes of 
the dimensionless driving force e, and a chaotic state at yet larger values of e. Our study 
focuses on the large aspect ratio limit in which the wavelength of the wave is much smaller 
than the lateral dimension of the fluid layer. 

Transitions to spatiotemporal chaos have been observed in a variety of physical systems, 
with recent efforts concentrating on the study of fluid systems [§]. The fact that physical 
properties of fluids are normally well characterized, and that in some cases the transition 
to chaos occurs at reasonably small values of the control parameter, makes fluid systems 
especially attractive for an accurate quantitative study of the chaotic state. Among these 
recent studies we quote rotating convection in a circular geometry , spiral defect chaos in 
Rayleigh-Benard convection [Q, and Faraday waves [^]. Theoretical work, on the other hand, 
has concentrated on the study of coupled lattice maps in which simple maps that are 

known to be chaotic are arranged on a lattice and their dynamics coupled. Complementary 
studies focus on phenomenological models such as the complex Ginzburg-Landau equation 
0. This equation models the low frequency behavior of a class of systems close to onset 
of instability, and exhibits various types of chaotic behavior. A third approach, which we 
follow here, concerns the so called order parameter equations. The best known example of 
the latter is the Swift-Hohenberg equation used to model Rayleigh-Benard convection in a 
Boussinesq fluid p. Its associated amplitude equation, the real Ginzburg-Landau equation, 
coincides with the amplitude equation directly derived from the equations governing fluid 
motion. The Swift-Hohenberg equation, supplemented with non-gradient terms, has been 
recently used to study spiral defect chaos in Rayleigh-Benard convection [p!0|]. 

Recent experiments on Faraday waves in large aspect ratio systems have revealed a 
number of interesting phenomena P,|TI1-[T5[] . Among them, periodic standing wave patterns 



near onset are found to be unstable against a TAM at some finite supercriticality P, [TT| , [l3 . 



Associated with this instabihty, temporal fluctuations of the pattern have been observed, 
with a characteristic time scale that decreases continuously with increasing e. Beyond a 
certain value of e, the wave patterns appear temporally chaotic and spatially disordered 
PJT^ . In the disordered regime, the amplitudes of a range of Fourier modes appear to exhibit 



Gaussian statistics (deviations from Gaussianity have also been recently reported [Q). 
Interestingly, disordered wave patterns in the chaotic regime have been found to exhibit 
highly ordered time averages [^. The origin of the TAM instability and its relationship 
with the disordered state was first studied theoretically by Ezerskii et al. ||TI1, who derived 
a one dimensional model to describe the modulation. Later, Milner [|17| derived a set of 
coupled amplitude equations for a two dimensional surface, including nonlinear interaction 
and damping of the waves. He showed that a pattern of square symmetry is realized near 
onset in fiuids of low viscosity, and that it can become unstable against a TAM. However, its 
finite threshold and its dependence on the detuning parameter are not yet well understood 
0. Studies of chaotic wave patterns were also conducted by Rabinovich and coworkers [ITB 



who numerically studied an amplitude equation for a pair of counter-propagating waves, and 
found features similar to the chaotic states observed in the experiments. 

The model that we have derived for a complex order parameter ip preserves the rota- 
tional symmetry of the fiuid, reduces to Milner's amplitude equations when if) is assumed 
to be modulated along particular directions of propagation, and to Ezerskii et al. one di- 
mensional equation in the case they considered. We confine ourselves in this paper to the 
results obtained when only the simplest nonlinear term is retained in the model. With 
this choice of nonlinearity we find a standing roll wave pattern close to onset. The exact 
form of the nonlinear term in momentum space, and the results corresponding to different 
choices of nonlinearity that yield standing waves of square symmetry at onset will be pre- 
sented elsewhere ||T^. We then obtain the stability boundaries of the roll pattern for either 
phase (Eckhaus and zig-zag) or amplitude (TAM) instabilities. A small but finite nonlinear 
damping coefficient appears to be essential for the existence of the finite driving threshold 



for the TAM instability. Unfortunately, and similarly to previous studies, we do not have a 
self-consistent procedure to introduce nonlinear damping into the model, probably reflecting 
the fact that at this level it has to be regarded as phenomenological. Numerical integration 
is finally used to study pattern selection above onset, and a transition to spatiotemporal 
chaotic states mediated by a TAM at finite e. 

The model equation that we study in this paper is, 

dti^ = -7V^ + ifr + j(i + v')v^ + 7n) iv-i'^, (1) 

where tlj{x,y,t) is a two-dimensional complex field. Its real and imaginary parts are linear 
combinations of the surface displacement away from planarity and the surface velocity po- 
tential of the irrotational part of the flow |T9| , pO| . Equation (|I]) has been made dimensionless 
by choosing 1/go as a unit of length, and 2/Q as a unit of time {Q is the angular frequency 
of the sinusoidal driving force, and go is critical wavenumber at onset, determined by both 
Q and the surface tension of the fluid). The linear damping coefficient is 7 = Auq^/Q, where 
u is the kinematic viscosity of the fluid. The nonlinear damping coefficient 7„ is of the order 
of 7, but its exact value is undetermined. The quantity e = {f — •y) / j is the distance away 
from threshold of the primary instability, where / is the dimensionless amplitude of the 
driving force. The phase e~*^*/^ in the motion of the surface has been scaled out of the fleld 
tlj{x,y,t). Equation (|l|) can be derived perturbatively from the fluid equations in the limit 



of weak viscous dissipation and small wave-steepness ||T9[. With the nonlinear term used in 
Eq. (|l]), standing roll wave patterns are observed just above 6 = 0. We believe, however, 
that our flndings about secondary instabilities, and the general features of the transition 
to spatiotemporal chaos are qualitatively similar to the more realistic case in which the 
standing wave pattern near onset exhibits square symmetry. 

For 6 > 0, the quiescent solution ip = loses stability. A new stationary roll solution 
can be found approximately by considering a one- mode Galerkin approximation iljQ{x,q) = 
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aqC*^" cos(gx) + 0{a^) with, 

^ q'-l- 477n/3 + - 1 - 477„/3)^ + (1 + 7^) [16(/^ - 7^)/9 - (g^ - 1)^] 

1+7^ 

sin 26*5 = (7 + 37„a;q/4)//, and cos2^^g = — 1 — a^)/ f. Another solution, which does 
not exist for q = 1 and is unstable with respect to uniform amplitude perturbations, will not 
be considered here. The neutral stability curve is eQ{q) = y'l + [3(1 — g^)/ (47)]^ — 1. The 
roll solution breaks the translational invariance of the base state. Hence a phase diffusion 
equation can be derived to study the stability of the roll solution. We find that the Eckhaus 
stability boundary is given by 
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- Inq' (27 + \lnq') i^-al + 1 - j , (2) 

and the zig-zag stable region is given by > |(g^ — 1). Fig. 1 shows the Eckhaus and zig-zag 
stability boundaries for 7 = 0.1 (this value of 7 is motivated by recent experiments [^ pT| , p!3[ ) 



and 7„ = 0.05. The reentrant shape of the Eckhaus boundary is a direct consequence of 
the existence of a small nonlinear damping coefficient 7„. In the small e limit {e <^ 7^), the 
region of stability is symmetric around g = 1 and is given by, e > -^^{^ ^ q^Yi whereas for 
7„ = 0, the stable region is e < ^^[1 + 3(1 — g^)/(87)]^ — 1, which lies entirely in the region 
g > 1. The parabolic stability boundary for £ <^ 7^ (7^ 7^ 0) can also be obtained from 
a standing wave approximation to Eq. (|I]), in analogy to Milner's calculation [|l^. This 
approximation, however, fails to reproduce Eq. (0) at larger values of e. 

Stability against transverse amplitude modulations can be studied by assuming that 
ip{x,y^t) = [1 + a{y,t)]il)Q{x,q)^ with a{y,t) = aoit) cos{Qy) {Q small) and linearizing the 
resulting equation for a{y,t). The TAM unstable region is given by, 

^(?^-l)^>^7X + 77n«J, (3) 

and is shown in Fig. 1 as the shaded region. Instability occurs for finite Q. If 7„ = 0, all roll 
solutions are unstable to TAM. If, on the other hand, 7„ 7^ 0, the region of wavenumbers 
around g = 1 is stable. At fixed 7, increasing 7^ increases the width of this stable region. 



We next turn to the results of our numerical calculation. We use a pseudospectral method 
to solve Eq. on a square grid of size GAtc x 6Att, with periodic boundary conditions. The 
number of Fourier modes used for each axis is 256. Time stepping is performed by a Crank- 
Nicholson scheme for the linear terms (including ip*), and a second order Adams-Bashforth 
scheme for the nonlinear terms. The time step used is At = 0.1. The initial condition, 
ip{x, t = 0) is a set of Gaussianly distributed random numbers, of zero mean and variance 
10"^. A typical run length is 10^ — 10^ time steps. 

Fourier modes with wavenumber close to g = 1 dominate in the early linear regime. As 
the system enters the nonlinear regime, the circularly averaged structure factor exhibits a 
dominant peak at g = Qmax which is seen to continuously shift away from 1 (detuning) and 
towards the Eckhaus stable region for all values of e, i.e., the dominant wavenumber of the 
roll pattern shifts to g > 1. The value of Qmax at long times is shown by the circles in Fig. 
(1). The error bars indicate the half-width at half-maximum of the peak. At e = ec(7n) 
the Eckhaus and TAM boundaries cross. For e < Ec, there is a region of stability against 
both Eckhaus and TAM perturbations. At sufficiently small e, an asymptotic stationary roll 
state is reached with q^ax inside this region. For sufficiently large values of e, qmax enters 
the TAM unstable region first. Therefore the reentrant shape of the Eckhaus boundary and 
the fact that it crosses the TAM boundary at finite values of e provides a mechanism for a 
finite threshold for the TAM instability. The shift in Qmax and the position of the Eckhaus 
and TAM boundaries depend strongly on the value of 7„. 

We finally describe the asymptotic temporal dependence of the configurations as a func- 
tion of e. The results reported correspond to 7 = 0.1 and 7^ = 0.05, but they are quali- 
tatively similar in other cases provided that 7n 7^ (for larger 7^, both the appearance of 
the TAM and the transition to chaotic states are seen at larger e). An almost perfect and 
stationary roll pattern is found for e = 0.02 (Fig. 2(a)). As e is increased to e = 0.05 a very 
slowly varying transverse modulation of the roll pattern is observed (Fig. 2(b)). The wave- 
length of the modulation is about 3.2Ao, with Aq = 2tt the critical wavelength at onset. The 
periodic and almost stationary modulation suggests that the finite amplitude modulated 
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pattern at e = 0.05 can be described by ipi^^ v) = ^o(l + bo cos{Qy)) cosqx + ■ ■ ■. Fig. 3(b) 
shows the zeros of ||'?/'|| {\\ip\\ < W'lpWmax/^ is considered a zero). They all originate from 
the cosqx factor. Hence the TAM is weak enough that ||1 + bocos{Qy)\\ > 0. At e = 0.1 
(Fig. 2(c)), the modulation becomes stronger and additional zeros of WipW appear (Fig. 3(c)). 
We call the regions in which the slowly varying component of WipW is zero a TAM defect. At 
this value of e the wave pattern becomes time dependent, with the temporal variation of the 
wave pattern coming mostly from the motion of TAM defects. At e = 0.3 (Fig. 2(d)), the 
length scale of the modulation becomes smaller, and the density of TAM defects increases 
(Fig. 3(d)). In order to characterize the temporal fluctuations of the wave patterns in this 
state, we have calculated the power spectrum of a time series of both ip a.t a fixed point in 
space and the amplitude of a fixed Fourier mode. Figure 4 shows the time series of i?e['^] 
at a fixed point (Fig. 4(a)), and of the Fourier mode g=0.75 (Fig. 4(b)), and their corre- 
sponding power spectra (Fig. 4(c) and (d)) for e = 0.50. The power spectra are broadband 
and decay as a power law of frequency with z ~ 4.0 for frequencies in the range of 
2 X 10^^ < 2uj/Q < 2 X 10^^. Similar power law decay of the power spectra is also observed 
for e = 0.3. The power law decay and the value of the exponent z are in good agreement 
with experimental results 
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FIGURES 

FIG. 1. Stability diagram for 7 = 0.1 and 7„ = 0.05 in the (g, e) plane: boundary of existence 

of stationary roll solutions (dashed line, S); q,nax, the peak of the circularly averaged structure 
factor calculated numerically (circles with error bars indicating the half-width at half-maximum of 
the peak); neutral stability curve (thin solid line, N); Eckhaus stability boundary (thick solid line, 
E); zig-zag stability boundary (long-dashed line, Z). The TAM unstable region is the shaded area. 

FIG. 2. Typical configurations of Re[ip\ at long times (shown in gray scale) for four values of 
e. The initial configuration was random. 

FIG. 3. Zeros of ||'^|| {\\ip\\ < W'lpWmax/^ is considered a zero) for four values of e at the same 
times as in Fig. 2. 

FIG. 4. (a) Time series of Re[^\ at a fixed point, (b) Time series of a Fourier mode 
Re i^{qx = 0.75, = 0) . (c) Power spectrum of the time series shown in (a), (d) Power spectrum 
of the time series shown in (b). The dashed line represents a power law with an exponent of -4. 
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